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1 . 


INTRODUCTION 


The Summer Research Program (SRP), sponsored by the Air Force Office of Scientific Research 
(AFOSR), offers paid opportunities for university faculty, graduate smdents, and high school students 
to conduct research in U.S. Air Force research laboratories nationwide during the summer. 

Introduced by AFOSR in 1978, this innovative program is based on the concept of teaming academic 
researchers with Air Force scientists in the same disciplines using laboratory facilities and equipment 
not often available at associates' institutions. 

The Summer Faculty Research Program (SFRP) is open annually to approximately 150 faculty 
members with at least two years of teaching and/or research experience m accredited U.S. colleges, 
universities, or technical institutions. SFRP associates must be either U.S. citizens or permanent 

residents. 

Hie Graduate Student Research Program (GSRP) is open annually to approximately 100 graduate 
students holding a bachelor's or a master's degree; GSRP associates must be U.S. citizens enrolled 
time at an accredited institution. 

The High School Apprentice Program (HSAP) annuaUy selects about 125 high school students located 
within a twenty mile commuting distance of participating Air Force laboratories. 

AFOSR also offers its research associates an opportunity, under the Summer Research Extension 
Program (SREP), to continue their AFOSR-sponsored research at their home institutions through the 
award of research grants. In 1994 the maximum amount of each grant was increased from $20,000 to 
$25,000, and the number of AFOSR-sponsored grants decreased from 75 to 60. A separate annual 

report is compiled on the SREP. 

The numbers of projected summer research participants in each of the three categories and SREP 
“grants” are usually increased through direct sponsorship by participating laboratories. 

AFOSR’s SRP has well served its objectives of building critical links between Air Force research 
laboratories and the academic community, opening avenues of communications and forging new 
research relationships between Air Force and academic technical experts in areas of national interest, 
and strengthening the nation's efforts to sustain careers in science and engineering. The success of the 
SRP can be gauged from its growth from inception (see Table 1) and from the favorable responses the 
1997 participants expressed in end-of-tour SRP evaluations (Appendix B). 

AFOSR contracts for administration of the SRP by civilian contractors. The contract was first 
awarded to Research & Development Laboratories (RDL) in September 1990. After completion of the 
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1990 contract, RDL (in 1993) won the recompetition for the basic year and four 1-year options. 


2. PARTICIPATION IN THE SUMMER RESEARCH PROGRAM 

The SRP began with faculty associates in 1979; graduate students were added in 1982 and high school 
students in 1986. The follow'ing table show s the number of associates in the program each year. 


YEAR 

SRP Participation, by Year 

TOTAL 


SFRP 

GSRP 

HSAP 


1979 

70 



70 

1980 

87 



87 

1981 

87 



87 

1982 

91 

17 


108 

1983 

101 

53 


154 

1984 

152 

84 


236 

1985 

154 

92 


246 

1986 

158 

100 

42 

300 

1987 

159 

101 

73 

333 

1988 

153 

107 

101 

361 

1989 

168 

102 

103 

373 

1990 

165 

121 

132 

418 

1991 

170 

142 

132 

444 

1992 

185 

121 

159 

464 

1993 

187 

117 

136 

440 

1994 

192 

117 

133 

442 

1995 

190 

115 

137 

442 

1996 

188 

109 

138 

435 ! 

1997 

148 

98 

140 

427 1 
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Beeinnin- in 1993 due to budget cuts, some of the laboratories weren’t able to afford to fund as many 
2E « in previous yean. Since .hen, the number of funded positions has remand fanly 

constant at a slightly lower level. 


3. RECRUITING AND SELECTION 

The SRP is conducted on a nationally advertised and compermve-selection basis .******* 
faculty and graduate students consisted primarily of the mailing of 8,000 52-page SRP brochures m 
chainLsonstf departments relevant to AFOSR research and to administrators of grants ^ accredited 
2” colleges, and technical institutions. Historically Black Colleges and Umversmes 
(THBCUsTand Minority Institutions (Mis) were included. Bmchures also went to all panictpanng 
USAF laboratories, the previous year's participants, and numerous indtvtdual requesters (over 1000 

annually). 

RDL placed advertisements in die following publications: Bfaci feues in of 

Change and IEEE Spearum. Because no participants list either Physics Today 
B&iering Nem as Wng their souree of learning about the program for the past several years, 
fSremfnts in these magazines were dropped, and die funds were used to cover mcreases tn 

brochure printing costs. 

High school applicants can participate only in laboratories located no more than 20 miles from their 
residence Tailored brochures on the HSAP were sent to the head counselots of 180 high schools in 
die vicinity of participating laboratories, with instmctions for pubtazmg the program ” 

High school students selected to serve at Wright Laboratory's Armament D‘re c, ora.e(^ltn Am Foret 
Bad, Florida) serve eleven weeks as opposed to die etgh. weeks normally worked by high school 
students at all other participating laboratories. 

Each SFRP or GSRP applicant is given a first, second, and third choice of laboratory. High school 
students who have more than one laboratory or directorate near their homes are also given first, 

second, and third choices. 

laboratories make their selections and prioritize their nominees. AFOSR then determines the number 
to be funded at each laboratory and approves laboratories' selections. 

Subsequently, laboratories use their own funds to sponsor additional candidates. Some selectees do 
not accept the appoinsnem, so alternate candidates are chosen. This mula-step P ^“,S 

results in some candxtaes being notified of their acceptance after scheduled deadlines Tbe tourl 
applicants and participants for 1997 are shown tn thts table. 


3 


1 1997 Applicants and Participants 

PARTICIPANT 

CATEGORY 

TOTAL 

applicants 

SELECTEES 

DECLINING 

SELECTEES 

SFRP 

490 

188 j 32 | 

(HBCU/MI) 

(0) 

(0) 

(0) 

GSRP 

202 

98 

9 

(HBCU/MI) 

(0) 

(0) 

(0) 

HSAP 

433 

140 

14 

TOTAL 

1125 

426 

55 


4. SITE VISITS 

During June and July of 1997, representatives of both AFOSR/NI and RDL visited each participating 
laboratory to provide briefings, answer questions, and resolve problems for both laboratory perscmnd 
and participants. The objective was to ensure that the SRP would be as constructive as possible for all 
participants. Both SRP participants and RDL representatives found these visits beneficial. At many of 
the laboratories, this was the only opportunity for all participants to meet at one time to share their 
experiences and exchange ideas. 


5. HISTORICALLY BLACK COLLEGES AND UNIVERSITIES AND MINORITY 
INSTITUTIONS (HBCU/MIs) 

Before 1993, an RDL program representative visited from seven to ten different HBCU/MIs annually 
to promote interest in the SRP among the faculty and graduate students. These efforts w ere marginally 
effective yielding a doubling of HBCI/MI applicants. In an effort to achieve AFOSR's goal of 10% 
of all applicants and selectees being HBCU/MI qualified, the RDL team decided to try other avenues 
of approach to increase the number of qualified applicants. Through the combined efforts of the 
AFOSR Program Office at Bolling AFB and RDL, two very active minority groups were found, 
HACU (Hispanic American Colleges and Universities) and AISES (American Indian Science and 
Engineering Society). RDL is in communication with representatives of each of these organizations on 
a monthly °basis to keep up with the their activities and special events. Both organizations have 
widely-distributed magazines/quarterlies in which RDL placed ads. 

Since 1994 the number of both SFRP and GSRP HBCU/MI applicants and participants has increased 
ten-fold, from about two dozen SFRP applicants and a half dozen selectees to over 100 applicants and 
two dozen selectees, and a half-dozen GSRP applicants and two or three selectees to 18 applicants and 
7 or 8 selectees. Since 1993, the SFRP had a two-fold applicant increase and a two-fold selectee 
increase. Since 1993, the GSRP had a three-fold applicant increase and a three to four-fold increase in 

selectees. 
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to addition to RDL's special recntiting efforts, AFOSR attempts each year to obtain additional funding 
or us^ 1 leftover funding cancellations the pas. year to fund HBCUM assonate, ^ year, 5 
HBCU/MI SFRPs declined after they were selected (and there was no one qualified to replace them 
with). The following table records HBCU/MI participation in this program. 



































1997 SRP FUNDING CATEGORY 

SFRP 

GSRP 

HSAP 

AFOSR Basic Allocation Funds 

141 

89 

_ 

USAF Laboratory Funds 

48 

9 

17 

-—-■—— 

HBCU/MI By AFOSR 
(Using Procured Addn'l Funds) 

0 

0 

N/A 

TOTAL 

9 

98 

140 


SFRP -188 were selected, but thirty two canceled too late to be replaced. 
GSRP - 98 were selected, but nine canceled too late to be replaced. 
HSAP - 140 were selected, but fourteen canceled too late to be replaced. 


7. COMPENSATION FOR PARTICIPANTS 

Compensation for SRP participants, per five-day work week, is shown in this table. 


1997 SRP Associate Compensation 


j PARTICIPANT CATEGORY 

1991 

1992 

r _ 

1993 

1994 

1995 

1996 

1997 

| Faculty Members 

$690 

$718 

$740 

$740 

S740 

$770 

$770 

Graduate Student 
(Master's Degree) 

$425 

$442 

$455 

$455 

S455 

$470 

$470 

Graduate Student 
(Bachelor's Degree) 

S365 

$380 

$391 

$391 

S391 

$400 

$400 

| High School Student 
| (First Year) 

$200 

$200 

$200 



$200 

$200 

High School Student 
(Subsequent Years) 

$240 



$240 

$240 

$240 

$240 


The program also offered associates whose homes were more than 50 miles from the laboratory an 
expense allowance (seven days per week) of $50/day for faculty and $40. day for graduate students. 
Transportation to the laboratory at the beginning of their tour and back to their home destinations at 
the end was also reimbursed for these participants. Of the combined SFRP and GSRP associates, 

65 % (194 out of 286) claimed travel reimbursements at an average round-trip cost of $776. 

Faculty members were encouraged to visit their laboratories before their summer tour began. All costs 
of these orientation visits were reimbursed. Forty-three percent (85 out of 188) of faculty associates 
took orientation trips at an average cost of $388. By contrast, in 1993, 58 % of SFRP associates took 
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orientation v isits at an average cost of $685 ; that was the highest percentage of associates opting to 
take an orientation trip since RDL has administered the SRP, and the highest average cost of an 
orientation trip. These 1993 numbers are included to show the flucmation which can occur in these 
numbers for planning purposes. 

Program participants submitted biweekly vouchers countersigned by their laboratory research focal 
point, and RDL issued paychecks so as to arrive in associates' hands two weeks later. 

This is the second year of using direct deposit for the SFRP and GSRP associates. The process went 
much more smoothly with respect to obtaining required information from tire associates, only 7% of 
the associates’ information needed clarification in order for direct deposit to properly function as 
opposed to 10% from last year. The remaining associates received their stipend and expense payments 

via checks sent in the US mail. 

HSAP program participants were considered actual RDL employees, and their respective state and 
federal income tax and Social Security were withheld from their paychecks. By the nature of their 
independent research, SFRP and GSRP program participants were considered to be consultants or 
independent contractors. As such, SFRP and GSRP associates were responsible for their own income 
taxes, Social Security, and insurance. 

8. CONTENTS OF THE 1997 REPORT 

The complete set of reports for the 1997 SRP includes this program management report (Volume 1) 
augmented by fifteen volumes of final research reports by the 1997 associates, as indicated below. 








APPENDIX A - PROGRAM STATISTICAL SUMMARY 


A. Colleges/Universities Represented 

Selected SFRP associates represented 169 different colleges, universities, and institutions, 
GSRP associates represented 95 different colleges, universities, and institutions. 

B. States Represented 

SFRP -Applicants came from 47 states plus Washington D.C. Selectees represent 44 states 
GSRP - Applicants came from 44 states. Selectees represent 32 states. 
jjSAP - Applicants came from thirteen states. Selectees represent nine states. 


| Total Number of Participants j 

SFRP 

189 

GSRP 

97 

HSAP 

140 

TOTAL 

426 | 
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SFRP Academic Titles \ 

Assistant Professor 

(A 

Associate Professor 


- 

Professor 

40 

Instructor 

0 


Chairman 


Visiting Professor 


Visiting Assoc. Prof. 


Research Associate 


TOTAL 


Source of Learning About the SRP 


Category 

Applicants 

Selectees 

Applied/participated in prior years 

28% 

34% 

Colleague familiar with SRP 

19% 

16% 

Brochure mailed to institution 

23% 

17% 

Contact with Air Force laboratory 

17% 

23% 

IEEE Spectrum 

2% 

1% 

BIIHE 

1% 

1% 

Other source 

KB 

8% 

TOTAL 

100% 

100% 
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APPENDIX B - SRP EVALUATION RESPONSES 


1. OVERVIEW 

Evaluations were completed and returned to RDL by four groups at the completion of the SRP. The 
number of respondents in each group is shown below. 


Table B-l. Total SRP Evaluations Received 


j- 

| Evaluation Group 

Responses 

j SFRP & GSRPs 

275 

| HSAPs 

113 

USAF Laboratory Focal Points 

84 

1 USAF Laboratory HSAP Mentors 

6 


All groups indicate unanimous enthusiasm for the SRP experience. 


The summarized recommendations for program improvement from both associates and laboratory 
personnel are listed below: 


A. Better preparation on the labs’ part prior to associates arrival (i.e., office space, 
computer assets, clearly defined scope of work). 

B. Faculty Associates suggest higher stipends for SFRP associates. 

C. Both HSAP Air Force laboratory mentors and associates would like the summer tour 
extended from the current 8 weeks to either 10 or 11 weeks; the groups state it takes 4- 
6 weeks just to get high school students up-to-speed on what’s going on at laboratory. 
(Note: this same argument was used to raise the faculty and graduate student 
participation time a few years ago.) 
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2. 1997 USAF LABORATORY FOCAL POINT (LFP) EVALUATION RESPONSES 


The summarized results listed below are from the 84 LFP evaluations received. 
I. LFP evaluations received and associate preferences: 



LFP Evaluation Summary. The summarized responses, by laboratory, are listed on the following 
page. LFPs were asked to rate the following questions on a scale from 1 (below average) to 5 (abo\e 

average). 

2. LFPs involved in SRP associate application evaluation process: 

a. Tune available for evaluation of applications: 

b. Adequacy of applications for selection process: 

3. Value of orientation trips: 

4. Length of research tour 

5 a. Benefits of associate's work to laboratory: 

b. Benefits of associate's work to Air Force: 

6. a. Enhancement of research qualifications for LFP and staff: 

b. Enhancement of research qualifications for SFRP associate: 

c. Enhancement of research qualifications for GSRP associate: 

7. a. Enhancement of knowledge for LFP and staff: 

b. Enhancement of knowledge for SFRP associate: 

c. Enhancement of knowledge for GSRP associate: 

8. Value of Air Force and university links: 

9. Potential for future collaboration: 

10. a. Your working relationship with SFRP: 
b. Your working relationship with GSRP: 

11. Expenditure of your time worthwhile: 

(Continued on next page) 
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12. Q ualit y of program literature for associate: 

13. a. Quality of RDL's communications with you: 

b. Quality of RDL's communications with associates: 

14. Overall assessment of SRP: 


Table B-3. Laboratory Focal Point Reponses to above questions 



AEDC 

AL 

L'SAFA 

PL 

RL 

WHMC 

TVL 

tt Evals Reev’d 


7 

1 

14 

5 

0 

46 

Question it 


2 

- 

86 % 

0 % 

88 % 

80 % 

- 

85 % 

2a 

- 

4.3 

n/a 

3.8 

4.0 

- 

3.6 

2b 

- 

4.0 

n/a 

3.9 

4.5 

- 

4.1 

3 

• 

4.5 

n/a 

4.3 

4.3 

- 

3.7 

4 

_ 

4.1 

4.0 

4.1 

4.2 

- 

3.9 

5a 

- 

4.3 

5.0 

4.3 

4.6 

- 

4.4 

5b 

- 

4.5 

n/a 

4.2 

4.6 

- 

4.3 

6a 

- 

4.5 

5.0 

4.0 

4.4 

- 

4.3 

6b 

- 

4.3 

n/a 

4.1 

5.0 

- 

4.4 

6c 

• 

3.7 

5.0 

3.5 

5.0 

- 

4.3 

7a 

- 

4.7 

5.0 

4.0 

4.4 

- 

4.3 

7b 

- 

4.3 

n/a 

4.2 

5.0 

- 

4.4 

7c 

- 

4.0 

5.0 

3.9 

5.0 

- 

4.3 

8 

- 

4.6 

4.0 

4.5 

4.6 

- 

4.3 

9 

- 

4.9 

5.0 

4.4 

4.8 

- 

4.2 

10a 

- 

5.0 

n/a 

4.6 

4.6 

- 

4.6 

10b 

- 

4.7 

5.0 

3.9 

5.0 

- 

4.4 

11 

- 

4.6 

5.0 

4.4 

4.8 

- 

4.4 

12 

- 

4.0 

4.0 

4.0 

4.2 

- 

3.8 

13a 

- 

3.2 

4.0 

3.5 

3.8 

- 

3.4 

13b 

- 

3.4 

4.0 

3.6 

4.5 

- 

3.6 

14 

- 

4.4 

5.0 

4.4 

4.8 

- 

4.4 
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3. 1997 SFRP & GSRP EVALUATION RESPONSES 


The summarized results listed below are from the 257 SFRP/GSRP evaluations received. 

Associates were asked to tale the following questions on a scale from 1 (below average) to 5 (above 
average) - by Air Force base results and over-all results of the 1997 evaluations are listed after the 

questions. 


1. The match between the laboratories research and your field. 

2. Your working relationship with your LFP: 

3. Enhancement of your academic qualifications: 

4. Enhancement of your research qualifications: 

5. Lab readiness for you: LFP, task, plan: 

6. Lab readiness for you: equipment, supplies, facilities: 

7. Lab resources: 

8. Lab research and administrative support: 

9. Adequacy of brochure and associate handbook: 

10. RDL communications with you: 

11. Overall payment procedures: 

12. Overall assessment of the SRP: 

13. a. Would you apply again? 

b. Will you continue this or related research? 

14. Was length of your tour satisfactory? 

15. Percentage of associates who experienced difficulties in finding housing. 

16. Where did you stay during your SRP tour? 

a. At Home: 

b. With Friend: 

c. On Local Economy: 

d. Base Quarters: 

17. Value of orientation visit: 

a. Essential: 

b. Convenient: 

c. Not Worth Cost: 

d. Not Used: 

SFRP and GSRP associate's responses are listed in tabular format on the following page. 
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Table B-4. 1997 SFRP & GSRP Associate Responses to SRP Evaluation 




Brooks 

Edwmf^i 


Griffis 

llansmai 1 


KirtLuxi 

LarkLiod 

Robins I 

TytxhiQ 1 

wpafb I 


1 

res 

H 

48 

4 


31 



33 

I 







mm 

ehi 

ES1I 

eh 

4.9 1 

4.6 

EH! 


EH 

EH 

EH 

HH 


gy 

eh 

EHI 

EOI 

EH 


IO! 

eh 

B3SS 

EH 

EH 

eh 

EH 


mm 

WWM 

eh 

EHI 

rwm 

4.2 1 

43 

EH 

IBiE 

EH 

EH 

EH 

EH 


mm 

ga 

wrm 

EH 

EH 

BE 

43 

EH 

KH 

EH 

EH 

EH 

EH 


mm 

WWM 

3 3 

EQ 

EH 

mm 

43 

MB 

K&E 

eh 


EH 

EH 

6 

4.3 1 

43 

k n 

eh 

EH 

KH 

4.0 

3.8 

— 

EH 


EH 

EH 

7 

■ 9 

gfl 

42 

4.8 

mm 

EH 

43 

4.1 

EHTE 

EH 

eh 

EH 

EH 

8 

111 

eh 

EU 

EO 

EH 


43 

go 


EH 

EH 

EH 

EH 

9 

mm 

in 

eh 

m 

EH 

4.5 

EO 

EH 

KH 



EH 

EH 

10 

mm 

fFl 

mm 

EO 

4.1 

4.1 

4.0 1 

EH 


EH 


EH 

43 I 

11 




gfil 

3.9 

4.1 

EO 

EH 


EH 

EH 

EH 

EON 

12 


mm 

4.3 

EO 

EH 

4.9 i 

EO 

4.6 

ESEI 

EH 

EH 

EH 

EHl 






Numbers below are 

pereentases 





13a 

83 

90 

83 

93 

87 

75 

100 

81 

100 

100 

100 

86 

87 

13b 

100 

89 

83 

too 

94 

98 

100 

94 

100 

100 

100 

94 

93 

14 

83 

96 

100 

90 

87 

80 

100 

92 

100 

100 

70 

« 

88 

15 

17 

6 


33 

20 

76 

33 

25 

HEdE 

100 

20 

8 

39 

16a 


26 

17 

9 

38 

23 

33 

4 

- 

- 

- 

30 


16b 

too 

33 

m 

40 

_ 

8 

- 

- 

- 

- 

36 

2 


16c 


41 

83 

40 

62 

69 

67 

96 

100 

100 

64 

68 


16d 



m 


• 

. 

- 

- 

- 

- 

- 

iHl 


17a 

, 

33 

100 

17 

50 

14 

67 

39 

- 

50 

40 

31 

35 

17b 


21 

m 

17 

10 

14 

. 

24 

- 

50 

20 

16 

16 

17c 

m 

m 

m 

• 

10 

7 

- 

- 

- 

- 

- 

2 

3 

17d 

100 

46 

- 

66 

30 

69 

33 

37 

100 

- 

40 

51 

46 
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4. 1997 USAF LABORATORY HSAP MENTOR EVALUATION RESPONSES 
Not enough evaluations received (5 total) from Mentors to do useful summary. 
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5. 1997 HSAP EVALUATION RESPONSES 


The summarized results listed below are from the 113 HSAP evaluations received. 

HSAP apprentices were asked to rate the following questions on a scale from 
1 (below average) to 5 (above average) 


1. Your influence on selection of topic/type of work. 

2. Working relationship with mentor, other lab scientists. 

3. Enhancement of your academic qualifications. 

4. Technically challenging work. 

5. Lab readiness for you: mentor, task, work plan, equipment. 

6. Influence on your career. 

7. Increased interest in math/science. 

8. Lab research & administrative support. 

9. Adequacy of RDL’s Apprentice Handbook and administrative materials. 

10. Responsiveness of RDL communications. 

11. Overall payment procedures. 

12. Overall assessment of SRP value to you. _ _ 

13. Would you apply again next year? Yes (9 c) 

14. Will you pursue future studies related to this research? Yes 

15. Was Tour length satisfactory? Yes (8 ' 


M 

umm 

■I 

El 

bhh 

esssbbsehii 

■Hj 

EZ5DBE3S1 


2.8 3.3 

4 4 4.6 

3.4 3.5 

4.5 4.8 

RH|| 

■ 

EflHBSS 

bi| 


1 

| 


mm 

Kpa 

m^Hi 

5 

6 

■ 

|ESHESi 

513a 


| 

mi 

7 

8 

B— 

—a 

BKES 

■ 

|| 

KHSI 

9 

in 

MEM— 

M 
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m 

MfMMm 

■ 

■ 

1 

1 

EEESI 
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Numbers below are percent! 

ises 


13 

14 

60% 95% 

70% 80% 

100% 100% 
71% 80% 

85% 100% 

54% 100% 

100% 100% 
71% 80% 

90% 92% 

65% 68% 

15 


ilia 

100% 50% 

msmmm 

80% 82% 
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INCORPORATING CONDENSATION INTO NASTD 


Jessica L. Thomas 
Graduate Research Assistant 
Department of Mathematics 
University of Tennessee Space Institute 

Abstract 

The goal of this study is to inexpensively obtain flow conditions which ideally are suited for flow visualization 
techniques involving vapor screens by incorporating condensation routines into the McDonnell Douglas NASTD 
code. The analysis follows that by Ryzhov, Pirumov, and Gorbunov and results in the Fokker-Planck equation 
which must be solved numerically by iteration. Further work will involve writing appropriate subroutines to 
perform this task and incorporating them into NASTD which will then be used to calculate the optimal flow 
conditions for flow visualization. 


y 
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INCORPORATING CONDENSATION INTO NASTD 
Jessica L. Thomas 

Introduction 

A popular method of flow visualization in wind tunnel testing involves vapor screens which rely on the presence of 
water droplets in the flow field of interest. Before running an expensive test the scientist would like to have a good 
idea of the initial conditions required to produce the desired droplet size and distribution during the actual test. To 
obtain these conditions quickly and inexpensively CFD codes with condensation capability can be used {1}. The 
overall purpose of this study is to improve the McDonnell Douglas NASTD code be incorporating routines to 
calculate condensation in the computational domain. To assist in the development a flow chart was made listing all 
subroutines called within the code. Then a literature search was performed to obtain a better idea of the physical laws 
governing condensation. 


Methodology and Results 

In this study the analysis follows that given in {2}. Interest is placed on spherical liquid droplets in a vapor medium 
with the assumption of equilibrium between the droplet and surrounding vapor. Also, the liquid and gaseous phases 
have the same temperature. To describe the nucleation of these spherical droplets, the change in the Gibbs free 
energy is used: 

AG = glu^p^T)-ii v (p v ,T)] + 4jrr 2 (j (1) 

where g is the number of molecules in a droplet with radius r, p./ and p v are the chemical potentials of the liquid 
and the vapor, p v and p^ are the vapor pressure and saturation vapor pressure, and a is the specific surface free 
energy (surface tension). The radius can be written as: 


y 




where V, is the volume per molecule. It is obvious from (2) that: 

g = 4Ttr i V l 


(3) 


To determine the equilibrium conditions, AG is minimized: 
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Aj) where the cluster becomes smaller. Also, the formation of larger clusters can result from the collision of two 
smaller clusters: 


A i + A j 


A i+j 


or can separate into two smaller clusters: 


A i+j 


A i+Aj 


However, except for the case of formation of dimers: 

Aj + Aj + Aj O A2 + Aj; Aj + Aj + B <=> A2 + B 

where B is a third substance such as the carrier gas, these reactions have little affect on condensation. Condensation 
occurs when the growth process dominates the decay process and the number of large clusters increases. The rate of 
formation of clusters containing g molecules Ig can be found if the reaction rate constants are known. The forward 
reaction constant is given by: 


K^P =a c 4nr 2 


Pv 


( 2nmkT) 


1/2 


( 10 ) 


where a c is the condensation coefficient. Since the reverse reaction constant is difficult to determine the ratio 


of the reaction rates is: 




K 


(f) 


n s - I n I K ( g r) 


( 11 ) 


where 

n g = n 1 exp\-—j ( 12 ) 

From the evaporation and condensation reactions, the number of clusters of class g per unit volume fg can be 
obtained by the following equation: 

dfg 

y -^- = i g -i g+ i,g = a . g) ( 13 ) 

where from (11) 
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I 


r _ t r (f)„ f _ v (r) f _ T,(f) fg-1 fg 

g K g n ifg-i K g h K g n i n g-i 

\ n g-I n g) 

For equilibrium conditions the number of clusters is given by ng and can be found from (12), but for the 
nonequilibrium conditions the system of differential equations (13) must be solved to find f CT . For very large g, 


where A g = ±1 (14) can be rewritten as: 


. -*•</>,„ ± f j 


/(s.t) = -K - n rt — - 
S 1 8 og n 


* D 1 °/ jAG 


Sg kT Bg 


where D= K n^ and Ag=dg. Another substitution of dl / og for I g - I g+l in equation (13) gives: 


= = d ( D ^) i 1 d \ D f dAG 

dt Bg dg dg) kT Bg 


This Fokker-Planck equation in future study is to be solved numerically for fg. 
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HUE ANALYSIS FACTORS FOR LIQUID CRYSTAL 
TRANSIENT HEAT TRANSFER MEASUREMENTS 


Derek E. Lang 
Graduate Student 
Department of Aeronautics 
University of Washington 

Abstract 

Hue Analysis is one approach for using liquid crystals to measure transient heat transfer 
in high speed flows. The indications given by wide-band crystals that are used to 
determine temperature, however, are highly sensitive to flow conditions, model material 
properties, viewing angles and illumination, and ffamegrabbing rate. This study used the 
hue capture range and duration of crystal colorplay as a measure the significance and 
influence of these factors in providing sufficient information to accurately measure 
temperature changes. An Experimental Design approach was used to define the test 
strategies and analyze the data. It was found that materials with high (pck) 1/2 , k/pc, and 
wide bandwidth crystals are the most critical factors. The hue analysis was then applied 
to various model geometries to illustrate the use of liquid crystals to measure heat 
transfer, and identify various boundary layer features. 
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HUE ANALYSIS FACTORS FOR LIQUID CRYSTAL 
TRANSIENT HEAT TRANSFER MEASUREMENTS 


Derek E. Lang 

Introduction 

The U.S. Air Force Academy/Aeronautical Research Center (ARC) is developing the 
capability to measure heat transfer in its Trisonic Wind Tunnel (TWT). One 
measurement technique assessed by ARC is liquid crystal thermography. This technique 
can be used to provide information about the global heat transfer on the surface of a the 
desired geometry. Thermochromic liquid crystals are surface coatings that change color 
with temperature. The monitoring of these temperature changes on a model surface with 
time are directly related to the heat transfer. Heat transfer measurements by liquid 
crystals are sensitive to a number of factors, ranging from the selection of the crystal and 
model materials to the application of the crystals and test conditions. This research 
analyzes the effect of these factors on liquid crystal measurements in the TWT test 
environments at ARC. 

Methodology 

The rate of heat transfer is dependent on the characteristics of the boundary layer formed 
at the model surface and the temperature differences driving the heat transfer. The heat 
flux equation below illustrates the dependence on the heat transfer coefficient, h, the 
adiabatic wall temperature, T aw , and the wall temperature, T w . 


q = h (T aw -T w ) 


The experimentalist determines heat transfer by measuring surface temperature(s) at 
given locations over time. Lang (Ref. 1) describes several methods for reducing this 
temperature data to heat transfer. 
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Ireland et. al. (Ref. 2) describe the liquid crystal measurement techniques used to 
accurately provide a single temperature measurement or continuous temperature data 
depending on the data reduction method used. In hue analysis, the hue of the liquid 
crystal is correlated with specific temperatures. Hue is one component along with 
saturation and value (or intensity) used to define liquid crystal color. The crystal 
temperature bandwidth is selected so that colorplay (i.e., hue changes in this case) occurs 
for the temperatures of interest. For example, a crystal bandwidth spanning the 
temperature range of the model surface during the test provides continuous inf ormation 
about the flow over the surface. 

Table 1 lists several factors contributing to the accuracy, or conversely uncertainty, of the 
heat transfer rate computed. This study focuses on those factors affecting the ability to 
measure temperature using liquid crystals. Experiments with liquid crystals on models in 
the TWT were conducted to evaluate effects of the models’ geometry and material 
properties (thermal conductivity, k, density, p, and specific heat, c p ), crystal bandwidth 
and total temperature. Model geometry effects are also related to heat-transfer rates, 
viewing angles, and illumination. 

Initially, the hues for the Hallcrest BM/R12C6W/C17-10 crystal with 6°C bandwidth and 
BM/12RC2W/17-10 crystal with 2°C bandwidth were correlated to temperatures in a 
static calibration rig. The Image Therm Engineering Liquid Crystal-Temperature 
Response Calibration Unit, MK II, insulated the liquid crystal which was sprayed onto a 
copper substrate and created a temperature gradient along the substrate length. Type-K 
thermocouples measured temperature at uniform distances along the substrate. The 
calibration rig was place inside the Trisonic Wind Tunnel to represent an in-situ 
environment. Florescent lights placed at 45-degree angles above and below the viewing 
window illuminated the substrate inside the tunnel. A Sony 3ccd Color Vision Camera 
(with RGB format) monitored the crystals at 30 Hz. A Matrox Meteor framegrabber 
captured the images at a user-defined rate and saved them on a personal computer. 
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Several models with a liquid crystal were then tested in the TWT with the test 
configurations in Table 2 to screen for significant factors under supersonic flow 
conditions. The TWT is a supersonic blowdown facility capable of producing up to 
Mach 4.28 flows. The TWT was operated at total temperatures of 22-35°C with a stilling 
chamber pressure of 1.2xl0 6 N/m 2 . The model surface temperature decreased to 
temperatures of -5°C once supersonic flow was established. Figure 1 shows the models 
used to represent the different heat transfer characteristics that ARC may investigate, 
including: sharp and blunted cones, a missile, a lifting body, and airframe-propulsion 
integrated vehicle forebody (i.e., API model). Table 3 lists the model materials and 
properties. Omega “Cement-On” type-E thermocouples were attached to the model 
surfaces to verify temperature readings. The liquid crystal was applied using airbrush 
with a minimum number of passes over the surfaces to ensure a thin coating. 

The measures of these factors’ influences were hue range captured by the video 
framegrabber as temperatures pass through the crystal bandwidth and the duration of this 
colorplay. Hue range is the usable portion of the hue spectrum that can be correlated to a 
temperature. Hue range also indicates the accuracy of the temperature measurement. 
Duration of colorplay relates to the amount of time available to capture the hues and the 
speed of ffamegrabbing required. The Experimental Design Method described by Smith 
and Launsby (Ref. 3) served as the basis for the test configurations selected and data 
analysis. 


Results 

Static Calibration. Figures 2a and 2b present the calibration curves for the 6°C and 2°C 
crystals. The lower end of the hue scale from 0-0.17 corresponds to the color red. The 
colors then change from yellow to green and then blue at the upper end of the scale. 
Initially, it was difficult to obtain data in the hue range between 0.15 - 0.40 which 
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encompasses the green, yellow, and red colors. As seen in the figures, the hues change 
rapidly with small increases in temperature for this portion of the hue spectrum. The 
calibration rig was limited in its ability to create the small temperature gradients 
necessary to isolate hues in this range. The colorplay from blue to red occurred within 
such a short distance along the length of the substrate that the ffamegrabber could not 
capture hues corresponding to the colors in between. Later, the crystals were cooled to 
temperatures below the crystal bandwidth at one end of the substrate and allowed to 
warm to ambient temperature. Colorplay was then spread across a wider portion of the 
substrate, and thermocouples recorded temperatures as the hue changes occurred. The 
temperature component of standard deviation from the best-fit curve was 0.55°C for the 
6°C crystal and 0.52°C for the 2°C crystal. Because the entire hue spectrum must be 
spanned within the crystal bandwidth, the narrower bandwidth crystals will have a flatter 
temperature-hue curve. As a result, the impact of the ±0.3°C thermocouple uncertainty of 
due to noise results in wider variances from the calibration curve-fit given uncertainty in 
hue determination is constant between the 6°C and 2°C crystal. An optimum bandwidth 
must be selected, since the wider bandwidth actually increases the uncertainty in the hue. 

The blue color is visible for several degrees above the manufacturer-specified upper 
temperature limit. As a result, the temperature-hue curve flattens above a hue of about 
0.60, and the uncertainty in corresponding temperature increases. Visually, blue and red 
lead to black on either side of the crystal bandwidth. However, hue is measured on a 
radial scale so that reds occur at hues approximately 0-0.15, but reappear for hues from 
about 0.80 to 1.0. High hue values (e.g., above .85) with low temperatures were given 
subtracted from 1.0, so that they corresponded with reds at the lower hue spectrum. 
However, it was difficult to define the correct hue cut-off (i.e., reds might occur with hues 
at 0.75 up to 1.00); and it is probably better to ignore hue readings above those 
corresponding to the green or green/blue portion of the hue scale. 


2-6 


Wind Tunnel Tests. The hues recorded by the framegrabber and associated temperatures 
measured by thermocouple during wind tunnel tests are compared with the static 
calibration curve in Figure 2a. The hues captured on the cone and one of the missiles are 
nearly constant, even though the wall temperature decreases with time; resulting in a hue 
range less than 0.1. In contrast, the lifting body matches the static calibration and spans 
the entire usable hue range, as it should if surface temperature passes completely through 
the crystal bandwidth. Although a large hue range does not necessarily imply an accurate 
reading, a small hue range may be an indication of inaccurate temperature measurement 
along some portion of the hues captured. 

Before assuming a problem with the hue capture, it must be realized that hue range may 
be limited because the surface temperature, as in some of the API model cases, never 
passed completely through the bandwidth, either because it reached adiabatic wall 
temperature or the heat transfer rate was so low. In most cases however, colorplay is 
completed by the end of the test run. 

Hue range is also related to the ability to capture data points within the crystal bandwidth 
which is related to colorplay duration. The high heating rates, material properties, and 
framegrabbing speed allowed the framegrabber to collected only a few points for the cone 
presented in Figure 2a, and resulted in a small hue range. Figure 3 shows the cone 
temperature passes through the crystal bandwidth much faster than the lifting body. This 
colorplay duration thus is the lower threshold for the framegrabber speed. Figure 4 is the 
result of Analysis of Variance (ANOVA) applied to the colorplay duration, and confirms 
that material, heat transfer, and total temperature are the major factors (note: bandwidth 
has some effect due to statistical association). The framegrabbing speed was found to be 
relatively significant to the hue range in Figure 5, and should be explored further. Shutter 
speed may affect the hue accuracy, but framegrabbing speed would seem to only affect 
the number of data points captured. 
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Figure 2a shows that even with an adequate framegrabbing speed, the hue measurements 
are subject other factors leading to inaccuracy. The small hue ranges generally coincided 
with tests that had faint or “washed out” colors which could not be picked up as well with 
the ffamegrabber as shown in Figure 6. Lighting and glare problems existed for all of the 
tests, but were not likely the cause since the models showed very distinctive colors when 
viewed while warming up in the test section. Another potential cause of washout may be 
the camera shutter speed which should be examined further. 

Figure 7 compares combinations of material (pck) 1/2 , bandwidth, and tank temperature 
(indicative of total temperature) for the Stycast lifting body and Plexiglas cones, which 
represent the wide span of the material properties in question (note: the (pck) 1/2 of 
Ertalyte is actually slightly smaller than Plexiglas). Bandwidth is the most influential. 
Material property and its combination with bandwidth are also significant; whereas tank 
temperature has little effect. From the prediction equation below, the hue range is 
maximized for large bandwidths and high (pck) 1/2 : 

Hue Range = 0.164 + 0.1(B) + 0.07(A) - 0.07(-A)(B) 

Where -1 < A < +1, for bandwidths between 2°C and 6°C 

-1 < B < +1, for (pck) 1/2 between Plexiglas and Stycast 

The analysis above assumed that given different model geometries were comparable, 
given similar heat transfer coefficients (h>100 W/m 2 in this case). Figure 5 compares the 
relative effects of k/pc for a constant cone geometry and crystal bandwidth. The 
optimum setting would be either a low k/pc with high heat transfer or high k/pc with low 
heat transfer. 

The color of the wide-band crystal is sensitive to viewing angle and illumination. In the 
tunnel, model viewing was complicated by glare and shadows on curved surfaces. In the 
case of the cones, the location of glare and shadows coincided with model curvature. A 
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comparison of Ertalyte models attempted to isolate the effects of the two factors as shown 
in Figure 8. The analysis suggested that the best combination for maximizing hue range 
was a shadowed area at an angle to the line of sight. This seems counter-intuitive, and is 
even more suspect since the hue range on the lifting body indicated an indifference to 
lighting (though the Stycast material has other favorable characteristics that may 
outweigh the lighting issue). It does appear that the hue ranges on the Ertalyte models 
was limited more by heating and material considerations than by lighting issues. 
Lighting angle may also be the cause of the “tongue”-shaped color pattern present on the 
lower side of the nose region of the conical models. The reduction of the images 
repeatedly captures hues in this pattern. Additional testing should be conducted to ensure 
that this is an illumination issue and not due to asymmetric flow. 

Flow Characterization. Once the appropriate selection of model geometries, materials, 
and test conditions are determined, liquid crystals provide substantial information about 
the flow over the model surface either through flow visualization or heat transfer 
distributions. For example, flow transition is generally indicated by a rise in heat transfer 
trend moving aft of the leading edge. This is readily seen by the liquid crystals on the 
API model where the beginning of transition is located at point of slowest color change in 
Figure 9, i.e., just before heat transfer increases. Figure 10 also illustrates the use of wide 
band liquid crystals to identify other flow features over the surface of the models, such as 
flow separation, transitional flow striations, and shock impingement. 

Heat transfer data reduction is more complex depending on the assumptions and 
techniques used. Cook-Felderman or Kendall-Dixon methods are well suited for the 
transient conditions of the TWT. However, the impacts of the start-up temperature spike 
during a run must be carefully considered. Moreover, use of these methods generally 
requires a minimum sampling frequency of 50 Hz, which was greater than the current 
video-framegrabbing capability. Babinsky (Ref. 4) supplements liquid crystal 
measurements with thermocouple data to obtain global surface heat flux by temperature 
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ratios. Figure 11 compares the heat transfer over the Missile model using this 
proportional heat transfer approach to application of complementary error function, and a 
narrow bandwidth, single-temperature method. 

A common problem in determining the analytical solution for heat transfer is the lack of 
knowledge of the state of the boundary layer (i.e., laminar, turbulent, or transitional) and 
boundary conditions. Separation regions and shock interactions make it difficult to 
predict recovery temperature, and the use of total temperature as a substitute in the heat 
transfer equations is complicated by the low enthalpies flow regime of the TWT where, 
total temperature is greater than wall temperature while adiabatic wall temperature lower 
than wall temperature. The advantage of the multiple, or continuous, data readings from 
liquid crystals is the ability to implicitly determine recovery temperature. 


Conclusions 

The wide-band liquid crystals provide global measurements on a time-continuous basis. 
Materials with high (pck) 1/2 and k/pc, lower heat transfer, wide crystal bandwidth, and 
high sampling frequency. Hue range is an indicator of the accuracy of the temperature 
measurement; while colorplay duration defines the minimum framegrabbing speeds 
required. Bandwidth and material selection are the most factors affecting these 
parameters. Crystal bandwidths should extend across the entire spectrum of wall 
temperatures anticipated during the test. 

The size of the hue range required is dependent on the data reduction technique. For 
example, only one data point is needed for the complementary error function method. 
For this approach, it is more critical to correlate the duration of colorplay to the video¬ 
framegrabbing rate. If a continuous data reduction technique is used (e.g., Babinsky), 
then hue range is related to the required temperature range and measurement uncertainty 
levels. For the tests conducted, adequate heat transfer rates were obtained in the cases 
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where the hue range was at a minimum 0.15-0.2 and colorplay lasted at least 10 seconds. 
Additional study is necessary to understand the influence of shutter speed, viewing angle 
and illumination. Once the model material and bandwidth are selected to give accurate 
temperature-hue relations, liquid crystal thermography provides significant data on the 
flow characteristics over complex configurations. 
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Table 1: Factors Affecting Heat Transfer Measurements in Trisonic Wind Tunnel 


Factor 

Measurement 

Technique 

Data 

Reduction 

Flow 

Physics 

Comment 

Crystal Bandwidth 

X 



Constrained by Available 

Crystals 

Crystal Temperature 

X 



Constrained by Available 

Crystals 

Material (rho.c.k) 

X 




Material (k/rho,c) 

X 




Model Geometry 

X 

X 

X 


Total Pressure 



X 


Total Temperature 

X 

X 

X 


Mach Number 

X 

X 

X 

Focus onM = 4.28 

Coating Thickness 

X 

X 


Kept relatively constant; assume 
negligible affect on heat transfer 

Illumination 

X 



Can keep constant once 
configured 

Viewing Angles 

X 




Start Time Uncertainty 


X 



Initial Surface Temperature 


X 



Data Reduction 

Assumptions 


X 


Assume h converges to constant 


Table 2: Test Configurations 


Geometry 

Material B 

andwidth (deg) 

lank lemp (P) | 

Sharp Cone 

Plexiglas 

2 


Sharp Cone 

Plexiglas 

2 

110 

Sharp Cone 

Plexiglas 

6 

90 

Sharp Cone 

Plexiglas 

6 

110 

Sharp Cone 

Ertalyte 

6 

110 

Sharp Cone 

Agloflon 

6 

110 

Blunt Cone 

Plexiglas 

6 

110 

Blunt Cone 

Ertalyte 

6 

110 

Blunt Cone 

Agloflon 

6 

90 

Blunt Cone 

Agloflon 

6 

110 

Missile 

Ertalyte 

6 

110 

Lifting Body 

Stycast 

2 

90 

Lifting Body 

Stycast 

2 

110 

Lifting Body 

Stycast 

6 

90 

Lifting Body 

Stycast 

6 

110 

API 

Plexiglas 

6 

110 

API 

Ertalyte 

6 

110 
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Figure 1: Models used (a) Sharp Cone; (b) Missile; (c) Slightly Blunted Cone; 
(d) Airframe Propulsion Integration (API); and (e) Lifting Body. 


Table 3 - Model Material Properties 


Material 

Density 

(kg/m 3 ) 

Specific Heat 
(J/Kg.K) 

Thermal 

Conductivity 

(J/Kg.K) 

sqrt(rho.c.k) 

(Ws ,/2 /m 2 .K) 

k/rho.c 

(l/m 2 s) 

Plexiglas 

1190 

1462 

.19 

575 

6.90e-8 

Ertalyte 

1201 

1075 

1.26 

501 

1.50e-7 

Agloflon 

2220 

981 

.34 

861 

1.56e-7 

Stycast 

2277 

904 

1.24 

1594 

6.03e-7 


2-13 







































Average Effect on Colorplay Duration (s 



Time (seconds) 

Figure 3: Cone Wall Temperature vs Time 



A C -AC B -BC ABC -AB 

Factor 


Figure 4a: Comparison of (pck) 2 , Bandwidth and Tank Temperature 
Effects on Colorplay Duration 
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Average Effect on Hue Rang Average Effect on Colorplay Duration (s 



Factor 

Figure 4b: Comparison of Heat Transfer, k/pc, and Framegrabbing Speed 
Effects on Colorplay Duration 



Hue 

Figure 5: Comparison of Heat Transfer, k/pc, and Framegrabbing Speed 

Effects on Hue Range 
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Average Effect on Hue Rang 




Figure 6: Colorplay Washout on Visual and Hue Images 


0.120 

0.100 

0.080 

0.060 

0.040 

0.020 

0.000 

B A -AB -AC C ABC -BC 



Factor 
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Average Effect on Hue Range 


Figure 7: Comparison of (pck) 1/2 , Bandwidth and 
Tank Temperature Effects on Hue Range 


0.160 
0.140 
0.120 
0.100 
0.080 
0.060 
0.040 
0.020 
0.000 

-AB B A 

Factor 

Figure 8: Comparison of Surface Angle land Illumination Effects 

on Hue Range 
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Figure 9: Identifying Boundary Layer Transition on API Model 










Figure 10: Liquid Crystal Indications of Boundary Layer Features 



Figure 11: Comparison of Data Reduction Methods For Missile Model 
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A SETUP FOR PHOTOASSOCIATION 
OF 

COLD, TRAPPED CESIUM ATOMS 


Bridget McGrath 
Department of Physics 
University of Colorado at Colorado Springs 


Abstract 


An existing Magneto-Optical Trap (MOT) setup was modified to allow for 
photoassociation studies. A photoassociation laser was directed towards the trapping cell, 
and a system for molecule detection was installed, consisting of a photoionization laser, 
and channeltron detection. Initial attempts at photoassociation were unsuccessful, but 
work will continue in this area, now that the groundwork has been completed. 
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Introduction 


In the past 10 years, successes in the field of laser trapping and cooling of atoms 
have furthered research into quantum effects and precision spectroscopy. As a next 
logical step, interest has turned to performing the same feat for molecules. This stems 
from the applications of cold molecules to developing better clocks and stabilized lasers. 
Laser cooling involved the repeated momentum exchange of many photons. However, the 
cycling transitions in the alkali atoms which allow for rapid absorption and emission of 
photons, necessary for cooling, are not available in molecules. Because of the molecules’ 
additional rotational and vibrational degrees of freedom, there are no transitions which can 
be easily exploited in this fashion. A complicated scheme has been proposed for possible 
cooling of molecules. [1] Most of the recent research, however, has focused on the 
creation of cold molecules by photoassociation of already cooled and trapped atoms. [2, 3, 
4, 5, 6] 

Photoassociation studies of Li, K, Na, and Rb have allowed spectroscopy of the 
long range molecular states largesly inaccesible from bound molecular states. In this paper 
we report the results of an attempt to create cold Cesium molecules in the same fashion.. 

The majority of photoassociation studies use atoms trapped in a Magneto-Optical 
Trap (MOT). A MOT is comprised of basically two features. The first, known as optical 
molasses involved 3 counterpropagating laser beams tuned slightly to the red of the 
desired transition in the atom. Doppler shifts arising from an atoms motion will bring the 
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laser light into resonance with the atom. Absorption of the photon gives the atom a 
momentum kick in the opposite direction of its initial velocity. The total effect is to slow 
the atoms within the laser beams. The second component of the MOT is a pair of 
Helmholtz coils. The magnetic field created by the current creates, via the Zeeman shift, a 
position-dependent force on the atoms. Thus the atoms are slowed and trapped. 

In the photoassociation process, a photon is absorbed by a cold atom in close 
proximity to another. This can form an excited state bound molecule which will decay 
90% of the time back into two atoms. Ten percent of the excited molecules will decay 
into a bound molecular ground state. Either case will result in a loss of atoms from the 
trap. Thus a photoassociation spectrum can be achieved by studying trap loss as a 
function of the frequency of the photoassociation laser, as has been in the previous studies. 
A second possible method for studying photoassotiation involves detecting the resultant 
molecules directly. 

These molecules can be ionized by the absorption of another photon, and the ions 
can then be detected using an electron multiplier. [7] Time of flight analysis could then 
provide the means of differentiating atoms from molecules. A photoassociation spectrum 
can also be achieved in this manner. A further benefit of this approach is that, since the 
molecules have many more intermediate states than the atoms, photoionization is much 
more likely process, and saturation can easily be achieved for relatively low ionization 
laser intensities. 
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Methodology 


Our setup consisted of a standard MOT with a background pressure of 10" 10 Torr. 
The trapping beams come from an 852 nm diode laser locked to the cycling 4-5’ transition 
by saturated absorption, with a series of acousto-optic modulators (AOM’s) providing the 
detuning of 20 MHz. The repumping light is provided by another diode laser on the 894 
nm F= 3-4 transition. Our magnetic field gradient is 10 Gauss/cm, and the MOT captures 
about 10 7 atoms at a density of 10 9 cm" 3 . 

In order to achieve a higher density, we modified the MOT by inserting dark spots 
into the center of the repumping beams, as first proposed by Ketterle et. al. . [8] This 
“dark spot “ MOT confines the atoms in the central region of the trap to the “dark” F=1 
hyperfine ground state, so that they cannot interact with the trapping laser, thereby 
reducing the outward radiation pressure at the center of the trap. 

A pulsed frequency doubled ND:YAG laser (532 nm) ionizes the atoms and 
molecules, and the subsequent ions are collected and detected with a Channeltron electron 
multiplier. 

Results 

We initially performed a series of tests on the setup. With a two- photon process 
required for the ionization of Cesium, we expected that the Channeltron signal strength 
should increase as E 2 of the ionization laser. As seen in Fig. 1, we found an experimental 
slope of 1.5, when plotting the log of the signal size vs. The log of the energy. 
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Fig. 1 

Ion Signal Size vs. intensity 
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Following that we focused on distinguishing the molecular signal from the atomic 
signal, in preparation for later attempts at photoassociation. We found that the Cs ions 
reached the detector between 25 and 30 microseconds after the laser pulse (Fig. 2)., in 
agreement with our computer simulation for the setup. Even without photoassociation, we 
expect a certain number (about 1%) of naturally occurring molecules in the trap. These 
would almost all be ionized at our laser power, and collected by the Channeltron. There is 
indeed a smaller peak occurring at a factor of 2 1/2 later in time of flight. Identifying this 
definitively as the molecular peak, however, gives rise to several problems. 
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Fig. 2 

Histogram of Pulse Frequency vs. Time of Flight 



Time of Flight (fis) 


First, there is also a third, much larger peak, occurring later, for which we have no 
solid explanation. It is possible that this is the result of ionization of Cesium from the 
walls of the cell, but we are far from certain. Secondly, as the intensity of the ionization 
laser is reduced, one would expect the atomic peak to fall off sharply, while the molecular 
peak would be relatively unaffected, due to saturation. This, however, was not the case. 
We found that when the intensity was reduced by a factor of 10, the frequency of 
ionizations per pulse for atoms decreased by a factor of 100, as predicted. The molecular 
peak, on the other hand, experienced a frequency reduction by a factor of 10,000. 

Finally, it was difficult to get consistent data. As mentioned above, the time of the 
atomic peak varied in different runs, making the expected time of a molecular peak 


3-7 




difficult to identify. We believe that part of the problem lies in the stability of our MOT. 
Further modifications to the setup, such as locking the 894 laser with saturated 
absorption, and redesigning the cell, should provide the needed stability. Attempts at 
photoassociation were unsuccessful. Nevertheless, the groundwork for further studies has 
been completed. 
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RELATIONSHIP BETWEEN GROWTH HORMONE 
AND MYELIN BASIC PROTEIN EXPRESSION IN VIVO 
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ABSTRACT 

Proper functioning of the mammalian nervous system requires myelination of neuronal axons. 
Myelination of the newborn mouse brain begins shortly after birth and is complete by about 20 days 
of age. Abnormalities of myelin basic protein (MBP) production have a direct impact on myelination 
as has been demonstrated in the shiverer mouse. Mice with less than 25% of normal MBP levels have 
aberrant myelination and brain development and demonstrate a characteristic tremor at 12 days of 
age. Growth hormone (GH) and insulin-like growth factor I (IGF-1) have been shown in vitro and 
in vivo to affect myelination. This project explores the effect of GH deficiency on myelin basic 
protein expression to determine whether the GH deficiency exacerbates MBP haploinsufficiency. 
These studies grew out of studies undertaken on children with 18q- syndrome who have only a single 
copy of the MBP gene, are hypomyelinated, and are also shown to suffer growth hormone deficiency 
or insufficiency. A hybrid mouse model which mimics these deficiencies of 18q- patients is being 
developed to explore the reltationship between GH and MBP in vivo. 
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RELATIONSHIP BETWEEN GROWTH HORMONE 
AND MYELIN BASIC PROTEIN EXPRESSION IN VIVO 

Donna M. Lehman 


INTRODUCTION 

Expression of the myelin basic protein (MBP) gene is critical for normal myelination of the central 
nervous system (CNS). This has been most clearly demonstrated by the neuropathology of the 
mutant shiverer (, shi/shi ) mouse in which a large portion of the MBP gene is deleted and MBP protein 
is essentially undetectable (1,2). Consequently, the CNS has virtually no compact myelin. These 
mutants exhibit a characteristic tremor (shiver) appearing at 12 days of age, coincident with the 
normal age at which myelination occurs in rodents, and die prematurely usually during status 
epilepticus. That MBP deficiency is responsible for the phenotype is supported by mutant rescue with 
a MBP transgene (3). 

Cross-beeding experiments utilizing mice transgenic for the MBP gene have established the level of 
MBP expression sufficient for normal myelinogenesis (4). Heterozygous shiverer mice produce about 
50% of the normal level of MBP mRNA and a proportionate amount of MBP protein. These animals 
exhibit a normal behavioral and morphological phenotype; however, they have minor biochemical 
changes in myelin composition. Mice which express less than 25% of MBP display the shivering of 
the dysmyelinating phenotype. These data indicate that the level of MBP expression influences the 
assembly of myelin by oligodendrocytes which suggests that factors which influence MBP expression 
also influence myelin synthesis and assembly. Additionally, factors which regulate oligodendroglial 
proliferation and differentiation may affect myelinogenesis. In this regard, there is much evidence, 
both in vivo and in vitro, that growth hormone (GH) and insulin-like growth factor 1 (IGF-1) enhance 
myelination (5,6,7,8). GH-deficient mice, such as the little (lit/lit) mutant which has an abnormality 
of the GH-releasing hormone receptor (Ghrhr) (9), are hypomyelinated; the hypomyelination is 
corrected by early postnatal GH treatment (10). 

Based on these observations, GH may be a potentiator of MBP production. To test this, we have 
begun to examine the impact of naturally occurring GH deficiency (lit/lit) on MBP in vivo, to 
investigate the effects of MBP deficiency ( shi/shi ) on GH production, and to explore the relationship 
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between GH and MBP in a hybrid mouse (lit/lit, +/shi) model. If GH affects MBP production, then 
MBP levels which are already reduced in the +/shi mouse may be further reduced in the GH-deficient 
hybrid, and therefore, the affected mice may display the dysmyelinating phenotype. Early treatment 
with GH may delay or prevent the appearance of the tremor. Alterations in developmental 
progression of myelination and its impact on function are being investigated. 

These studies in mice grew out of studies undertaken on children with 18q- syndrome (11,12,13). 
The characteristic features of this syndrome are short stature, mental retardation and deafness. 
Affected children have only a single copy of the MBP gene and are hypomyelinated as demonstrated 
by quantitative magnetic resonance relaxometry (14). Half of affected children are GH-deficient and 
the other children have evidence of dysregulation of GH production. Accumulating human data 
suggest that GH treatment may improve both myelination and cognitive function (unpublished) raising 
the possibility that a relationship exists between MBP and GH. These mouse studies permit the 
opportunity to more thoroughly explore the relationship between GH and MBP within the context 
of the developing brain and may provide insight into novel treatment approaches to children with 
hypomyelination. 

METHODOLOGY 

Mouse colony. Little and shiverer mutant mice were supplied by Jackson Laboratory. The 
little mutant is in the C57B1/6J background, and the shiverer mutant is in the C3H background. Wild- 
type little males were crossed with lit/lit females, and the heterozygous offspring (FI) were interbred 
to create +/+, lit/+, and lit/lil genotypes. C3H shiverer mice were also interbred to produce C3H 
mice for study. The hybrid mouse {lit/lit;shi/+) was produced by first crossing C57B1 lit/lit females 
with C3H shi/shi males to produce a double heterozygote. Male FI mice were backcrossed to C57B1 
lit/lit females. 

Mutant genotyping. A method was developed to genotype the little locus (ghrhr) by PCR and 
subsequent restriction analysis (15). The shiverer locus (MBP) was determined by Southern analysis 
as previously published (16). A 581-bp probe which spanned the deletion breakpoint was designed 
for this purpose. 
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MBP quantitation. MBP mRNA and protein levels were assayed in each mouse. Brains were 
removed from mice; the right hemisphere was processed for RNA and the left hemisphere was frozen 
in liquid nitrogen and stored at -70 °C for protein analysis by RIA. MBP mRNA levels were 
determined by Northern analysis. Total RNA was isolated from brain tissue using a protocol adapted 
from Chomczynski and Sacchi (17). Five pg of total RNA were electrophoresed through a 1.2% 
agarose gel. The blot was hybridized with a radiolabeled cDNA probe for the 14kD form of mouse 
MBP (gift of Dr. Carol Readhead). An a-tubulin probe was subsequently used as an internal control 
for normalizing the amount of RNA loaded. Quantitation was done with a phosphorimager. 

IGF-1 and GH quantitation. Blood was collected from mice by cardiac puncture. Serum was 
separated from the cellular components and frozen. Fifty pi of serum was used per assay; assays were 
performed in duplicate. Each sample was acid-extracted using C-18 cartridges (Sep-Pak, Waters 
Associates) and IGF-1 assayed by radioimmunoassay. Serum assay of IGF-1 was performed by Dr. 
Ross Clark at Genentech, Inc. Anterior pituitaries were harvested for analysis of intrapituitary GH 
and stored at -70 °C in sterile saline. 

Auditory Brainstem Evoked Responses (ABR). ABRs were recorded using a Bio-Logic 
Navigator system. The mice were anesthetized with inhaled 1.5% isoflurane in 60% 0 2 . One 
hundred microsecond clicks are generated at a repetitive rate of 11,4/sec as the sound stimulation and 
then delivered by earphone to an external ear canal of the mouse through a 20 cm connecting tube. 
The evoked potentials, taken with needle electrodes inserted at the vertex and the ipsilateral paw, 
were amplified X 10 4 , filtered with a band pass between 150 and 1500 Hz, and averaged from a total 
of 1000 click-evoked responses for the first 10 msec period following stimulation. The minimum 
detectable values on the two different points of the latency, measured by the cursor mode of the 
system, is 0.04 msec. The minimal hearing thresholds can be determined by 5 dB step sound 
stimulations according to the visually detectable amplitude of the first peak. 

RESULTS 

MBP levels in brain from little mice 

Although there is extensive data on the effect of GH/IGF-1 on myelination, there is no information 
on the specific effect of GH/IGF-1 on MBP in vivo. The period of maximal MBP expression in mice 
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is between 10 and 30 days after birth with the peak at about 15 days (18). Thereafter, express*, 
declines and plateaus a. a much lower level. Therefore, MBP expression was determined at 18, 25 
and 60 days postnatal in little mice. The data for 25-day mice is complete and show that ht/ht nuce 
have less MBP mRNA than controls (p < 0.05, + A vs. lit/lil) (Fig 1). This result suggests that a 
mechanism by which GH affects myelination is by altering the amount of MBP mRNA. The manner 
in which this is accomplished is unknown. MBP protein assays are in progress. 



Figure 1. Quantitation of MBP mRNAin 25-day lit 
mice. Data are shown as a ratio of MBP mRNA to 
ff-tiihulin mRNA. 


F^inrtional evaluation o f mvelination 

Alterations in the content or quality of myelin should have an impact on neural fimction. Auditory 

brainstem evoked responses were analyzed in several 2-month old mice at 90, 80, 70 60, 50, 45 and 

40 dbl: ABRs for 2 controls (C57B1/6 +/+), 2 heterozygotes (+/IU) and 2 littles {ht/ht) at 90 dbl are 
shown in figure 2. Wave responses were also recorded from homozygous shiverer ( shi/sht ) mice for 
comparison. These results show that there are distinct differences between the controls, 
heterozygotes and little animals. Analysis of the evoked responses revealed a Wave 1 in Ht/ht nuce 
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that had a normal amplitude, latency and morphology; in all animals tested, the remainder of the 
waveform lacked identifiable peaks. These data suggest that there is normal cochlear and VTIIth 
nerve function and abnormal auditory brainstem function in lit/lit mice. Whether the differences in 
ABR are due to hypomyelination is not yet clear; however, it is intriguing to consider that there is a 
relationship between GH deficiency, hypomyelination and hearing. It is noteworthy that homozygous 
shiverer and homozygous little mice also have distinct responses. This could indicate that either their 
degrees or basis of hypomyelination differ. Ultrastructural assay of the auditory nerve, such as by 
electron microscopy, may provide answers. Nonetheless, this functional assay can distinguish these 
genotypes which is helpful when studying the hybrid mouse model {lit/lit;shi/+). 


Testing at 70 dB nHI 



Figure 2. Auditory brainstem evoked responses in lit 
mice. 


Circulating IGF-1 levels in shiverer mice 

No information is available on the impact of the shiverer mutation on the neuroendocrine axis. 
Because of interest in the potential effect of an abnormality of MBP on GH production, serum IGF-1 
of 42-day old shiverer and control mice was assayed (Table 1). Due to the pulsatile nature of GH 
secretion, random serum GH determination does not provide useful information; therefore, IGF-1 was 
chosen as an initial marker for GH production. Since the number of 42-day old control animals is 
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small, statistical evaluation could not be undertaken; however, a trend is apparent. Assays on 
additional animals may make these data significant. The functional significance of these changes on 
parameters other than somatic growth has not been evaluated. Pituitaries have been gathered for 
analysis of intrapituitary GH. 


Genotype 

shi/shi 

+/shi 

+/+ 

IGF-1 

(ng/ml) 

421±104 

552 ± 135 

650 ± 99 

(n) 

(10) 

(13) 

(2) 


Table 1. IGF-1 levels in shiverer mice and controls 


Development of hybrid mouse model 

We have demonstrated that the hybrid {lit/lit; +/shi ) is viable. Studies are being conducted to 
establish the level of MBP expression and myelination in these hybrids and to determine if this differs 
from the +/+; +/shi and the lit/lit; +/+ mice during development. The functional consequences of 
these changes will also be investigated. 

DISCUSSION 

GH-deficient little mice have previously been shown to be hypomyelinated (10), but the mechanism 
remains unclear. Our data indicate that MBP deficiency may be a factor. Homozygous little mice 
have only about 72% of MBP mRNA levels of control mice. MBP protein assays are in progress. 

Auditory evoked brainstem responses have been used as a functional measure of myelination in 
shiverer mutants (19) as well as humans. We’ve shown that ABRs of little mutants are abnormal in 
that signal processing beyond the cochlea is disrupted. Very little data is available concerning GH- 
deficiency and auditory neurosensory impairment. This data may provide insight into the cause of 
abnormal ABRs in some 18q- children as well as potential therapy. Ultrastructural imaging of the 
myelinaton of the auditory nerve in little mice should prove useful in determining the cause of the 
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aberrant ABRs 


18q- patients are haploinsufficient for MBP, but their cause of GH-insufficiency is unknown. The 
GH-insufficiency region on chromosome 18 has been narrowed (12), and the MBP gene lies within 
this region. To determine whether MBP could be a candidate gene for GH-deficiency/insiifficiency 
in 18q- syndrome, we measured serum IGF-1 and intrapituitary GH levels in shiverer mutant mice 
and controls. The IGF-1 data indicate that homozygous shiverer mice have reduced circulating levels 
of this factor at 42 days. Since the number of control animals is small, statistical significance could 
not be determined. Assays on additional animals are underway. Also, intrapituitary GH analyses are 
in progress. 

Development of a hybrid mouse model which has only one copy of the MBP gene in a GH-deficient 
background should prove useful in examining the relationship between GH and MBP production in 
vivo. We have shown that this hybrid is viable and are initiating studies to determine the 
consequences of GH-deficiency on MBP haploinsufficiency. 
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